Resonant state expansion applied to two-dimensional open optical systems 
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The resonant state expansion (RSE), a rigorous perturbative method in electrodynamics, is ap- 
plied to two-dimensional open optical systems. The analytically solvable homogeneous dielectric 
cylinder is used as unperturbed system, and its Green's function is shown to contain a cut in the 
complex frequency plane, which is included in the RSE basis. The complex eigenfrequencies of modes 
are calculated using the RSE for a selection of perturbations which mix unperturbed modes of dif- 
ferent orbital momentum, such as half-cylinder, thin-film and thin-wire perturbation, demonstrating 
the accuracy and convergency of the method. The resonant states for the thin-wire perturbation 
are shown to reproduce an approximative analytical solution. 

PACS numbers: 03.50.De, 42.25.-p, 03.65.Nk 



I. INTRODUCTION 

The electro-magnetic spectrum of an open optical 
system is characterized by its resonances. Optical 
resonatorSfi such as planar microcavities^ photonic crys- 
tal fiberS)^ dielectric micro-spheres?^ micro-disksSiS and 
micro-cylinders^^— are examples of such systems designed 
to have a series of narrow resonances in their optical spec- 
tra. Such resonances, known as cavity modes in planar 
microcavities and whispering gallery modes (WGMs) in 
micro-spheres and -cylinders, are characterized by their 
spectral positions and linewidths, given, respectively, by 
the real and imaginary part of the complex eigenfrequen- 
cies of the system. 

Finite linewidths of resonances are typical for open 
systems and are due to energy leakage from the system 
to the outside. This leakage can be enhanced by var- 
ious structural imperfections and scatterers. In partic- 
ular, when an object is placed inside or in close prox- 
imity to the cavity, the resulting modification of the 
electromagnetic susceptibility perturbs the cavity reso- 
nances, changing both their position and linewidth, most 
noticeably for the high-quality (i.e. narrow-line) res- 
onances. This effect is the basis for resonant optical 
biosensingi^iiSiii The changes in the spectral properties 
of resonators in the presence of perturbations can be 
used to characterize the size and shape of the attached 
nanoparticlesii^ii^ The WGM resonances in microdisks 
and spherical microcavities have been used in sensors for 
the characterization of nanolayersi^., DNA^^ and protein 
molecules as well as for single atomic and nanoparticle 
detect ion F^^ii^ Other applications of high-quality modes 
include miniature laser sources^Sri^ and photonic-crystal 
optical fibers.^ 

While the eigenmodes of resonators with simple and 
highly symmetric geometries can in some cases be calcu- 
lated exactly, determining their perturbations presents a 
significant challenge as the popular computational tech- 
niques in electrodynamics, such as the finite difference 
in time domain^ii^ or finite element method^^"— need 
excessively large computational resources both in mem- 
ory and processor usagej ^^'^° This is due to the ex- 



tremely large computational domain in space and time 
required to model high quality modes. To treat such 
narrow resonances and their perturbations, we have re- 
cently developed^i a rigorous perturbation theory called 
resonant state expansion (RSE) and applied it to pla- 
nar optical resonators with different perturbationsH as 
well as to spherical resonators reducible to effective one- 
dimensional (ID) systems.— We have demonstrated on 
exactly solvable examples in ID that the RSE is a reliable 
tool for calculation of wave numbers and electro-magnetic 
fields of resonant states (RSs),— as well as transmission 
and scattering properties of open optical systems. 

In the present work we report an important step in 
the development of the RSE, extending the method to ef- 
fectively two-dimensional (2D) systems (i.e. 3D systems 
translational invariant in one direction) which are not re- 
ducible to effective ID systems. We provide a summary 
of the theory in a general 3D case, apply it to an effective 
2D system, and then consider examples which cover sev- 
eral important types of perturbations. We treat an ideal 
dielectric cylinder with uniform dielectric constant in vac- 
uum as unperturbed system and calculate perturbed RSs 
for homogeneous (i.e. reducible to ID) and inhomoge- 
neous perturbations, including a half-cylinder, thin-film, 
and thin-wire perturbations. None of these inhomoge- 
neous perturbations have known exact analytic solutions 
which could be used for verification of the RSE. How- 
ever, the case of a narrow wire inside a cylinder allows 
for an approximate analytic solution suitable for weak 
perturbations. ■^■^■■^^ This solution is compared with the 
present results of the RSE, demonstrating a good agree- 
ment. 

In the literature, the scattering properties of an open 
system are often described in terms of a continuum of its 
eigenstates, all having real frequencies. Such a contin- 
uum, together with isolated eigenmodes (e.g. waveguide 
modes), if they exist in the spectrum, form a mathemat- 
ically complete set of states suitable for expansions. The 
continuum, however, presents a significant obstacle for 
computational methods like perturbation theory, which 
are based on such expansions. The concept of RSs, natu- 
rally following from the observation of resonances in the 
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spectra of open systems, introduces another complete set 
of eigenstates which ehminates the continuum from the 
spectrum, replacing it by a countable number of discrete 
modes with complex frequencies. This is the basis used 
in the RSE. 

One striking feature of 2D systems, which is revealed 
in the present work, is the presence of a one-dimensional 
continuum in the manifold of RSs. This continuum is 
specific to 2D systems and is required for the complete- 
ness of the basis and thus for the accuracy of the RSE 
applied in 2D as discussed in this work in some de- 
tail. Such continua are generally known in the theory 
of quasi-guided modes in photonic crystal structures^ 
as potential sources of Wood-Rayleigh anomalies in opti- 
cal spectra.—"— Technically, they are caused in that case 
by the presence of square roots in the photon dispersion 
of light propagated and Bragg scattered inside the pho- 
tonic crystal. In the case of a dielectric cylinder, which is 
used for the RSE in 2D, the continuum originates mathe- 
matically from the cut in the cylindrical Hankel functions 
solving Maxwell's equations outside the cylinder. 

The paper is organized as follows. In Sec. II we give the 
general formulation of the RSE for an arbitrary three- 
dimensional (3D) system. In Sec. Ill we treat the special 
case of effective 2D systems using the homogeneous di- 
electric cylinder as unperturbed system and adding the 
contribution of the cut to the RSE. This is followed by 
examples illustrating the method and comparing results 
with existing analytic solutions. Details of the general 
formulation of the method, its application in 2D, and the 
calculation of matrix elements for specific perturbations 
are given in Appendices A-D. 



II. RESONANT STATE EXPANSION 



we find that E(r) = E,*(r) and fc = ±fc* also satisfy 
the same equation. Only — fc* has a negative imaginary 
part as required for time-decaying solutions. We label 
the resulting RS with index — n, so that fc_„ = 

For any open system, the RSs form an orthonormal 
complete set of eigenmodes. It follows from Eq. ([2]) that 
solutions decaying in time grow exponentially in space 
as r — >■ cx). Therefore the normalization of RSs cannot 
be simply given by the usual volume integrals over their 
wave functions but it also needs to involve the electro- 
magnetic energy flux through a surface surrounding the 
system. The orthogonality of RSs for n ^ m has the 
form 

= (kl-kl) f dre(r)E„(r).E„,(r) 
Jv 



dS E„ 



~d7 



E, 



ds 



(3) 



where the first integral in Eq. ([3|) is taken over an arbi- 
trary volume V which include all system inhomogeneities 
of £(r) while the second integral is taken over the surface 
Sy surrounding the volume V and contains the gradients 
d/ds normal to this surface. This follows strictly from 
Maxwell's equation Eq. ([T|) and Green's theorem, using 
V • E„ = on the surface Sv which is situated outside of 
the system inhomogeneities. It is convenient, following 
Eq. ([3]), to normalize the RSs in the way 
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In general, RSs in an open optical system with local 
dielectric constant £(r) and permeability = 1, where r 
is the three-dimensional spatial position, are the eigen- 
solutions of the Maxwell wave equation. 



Vx VxE„(r) = fc2£(r)E„(r), 



(1) 



which satisfy the outgoing wave boundary condition 

E„(r) ^ ,.-(I5-l)/2g»fc„r fQj. ^ ^ ^ (2) 

where r — |r|, D is the space dimensionality, fc„ is the 
wave-vector eigenvalue of the RS numbered by the in- 
dex n, and E„(r) is its electric field eigenfunction. The 
time-dependent part of the RS wave function is given by 
exp(— iaj„t) with the complex eigenfrequency ujn = ck„, 
where c is the speed of light in vacuum. RSs are ei- 
ther stationary or time-decaying solutions of Maxwell's 
equation. The wave numbers fc„ of time-decaying RSs 
lie in the lower half of the complex fc-plane and come in 
pairs, having the opposite real and equal imaginary parts. 
Indeed, if E„(r) and fc„ corresponding to RS n satisfy 
Eqs. (HI) and taking the complex conjugate of Eq. ([1]) 



where we have used in the second integral an ana- 
lytic continuation E{k, r) of the RS wave function E„(r) 
around the point k„ in the complex fc-plane. For solu- 
tions decaying in time both surface and volume integrals 
in Eqs. ^ and dU diverge for F — >■ oo. Their super- 
position, however, removes the divergencies making the 
normalization independent of V. In the special case of 
stationary states which are decaying in space, the surface 
integrals vanish for V ^ oo and the usual orthonormality 
in terms of infinite- volume integrals is restored. 

The RSE is based on the following three key elements. 
The first one is the Dyson equation, 

gfe(r,r') - G,(r,r') (5) 
-e I Gfc(r,r")Ae(r")afc(r",r')dr", 



which relates the perturbed and unperturbed GFs, Qk 
and Gfc, respectively. The difference between the per- 
turbed and unperturbed systems is a perturbation of the 
dielectric constant A£(r) with compact support. 

The second key element is the spectral representation 
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of the GF, 



Gfe(r,r') = 



;r)(g)E„(r') fc„ 



k(k — k„) w. 



(6) 



which takes into account simple poles of the GF aXk — kn 
as a sum and a cut of the GF in the complex fc-plane 
as an integral. Details of the derivation of Eq. ([6]) in a 
general 3D case, using the Mittag-Leffler and reciprocity 
theorems, are given in Appendix \^ accounting for the 
tensor form of the GF of Maxwell's equation. The pres- 
ence of the cut is a specific property of 2D systems having 
a continuum of RSs in their spectrum. The nature of the 
cut and its contribution to the RSE is discussed in more 
detail in Sec. llIII and in Appendix IB] The spectral repre- 
sentation Eq. ^ is used in the Dyson equation ^ , for 
both unperturbed and perturbed systems, equating the 
residues at the perturbed poles of both sides of Eq. ([5]) . 
The spectral representation of the unperturbed GF has 
the normalization constants Wn = '2kn following from the 
normalization condition Eq. as shown in AppenciesEl 
and[Bj As for the spectral representation of the perturbed 
GF Qk, in which E„ and fc„ are replaced by £i, and >r,y, 
it is not required by the RSE that the perturbed RSs are 
normalized in the same way, and thus the corresponding 
normalization constants can be any. 

The third element is the completeness of RSs, which is 
discussed in Appendix \X\ and mathematically expressed 
by the closure relation Eq. (jA9[) . It allows to expand the 
perturbed wave functions into the unperturbed ones: 



(7) 



The expansion coefficients Cm, are scaled here by 1/ \/k^, 
in order that the Dyson equation reduces to a linear sym- 
metric matrix eigenvalue problem^ 



(8) 



in which the matrix elements of the perturbation are de- 
fined by 



Vn,n = J Ae{r) E„(r) • E„.(r) dr . 



(9) 



Equations (I7]|ni), together with the normalization condi- 
tion Eq. Q used for the unperturbed modes, define the 
method called the RSE. 



III. APPLICATION TO 2D SYSTEMS 

Let us consider systems in 3D space which are homo- 
geneous in one direction (along the unit vector z of the 
z-axis), and thus can be reduced to effective 2D systems, 
as their wavevector component along z is conserved and 
the solution can be separated into a plane wave exp{ikzz) 



and the remaining (x, y)-problem which we express be- 
low in polar coordinates p — (p, tp) . For such a system, 
the solutions of Maxwell's equations split into two groups 
with orthogonal polarizations, called transverse electric 
(TE) and transverse magnetic (TM), where TE (TM) 
states have a electric (magnetic) field orthogonal to z. 
This nomenclature relates to the theory of waveguides 
where the light propagating along z has the dominant 
component kz of the wave vector, and the electric (mag- 
netic) field in TE (TM) modes is thus approximately per- 
pendicular to the wave vector. Although we restrict our 
treatment here to the opposite limit of kz = 0, we follow 
these adopted notations. In our case, however, the TE 
(TM) states have the magnetic (electric) field polariza- 
tion vector strictly parallel to z and thus normal to the 
wave vector of light at large distances. 

We treat in this work only the kz = TM states for 
which E„ = zEn and Eq. ([Ij reduces to 



ld_ 
pdp 



1 

p2 difi^ 



■e{p,ip)kl 



(10) 

The states are normalized according to Eq. Q in which 
the volume V is given by an infinitely long cylinder of 
radius R. In order to make the normalization constants 
finite, RSs are normalized per unit length along z. The 
normalization following from Eq. Q is given more explic- 
itly by^ 



1 



dip / pdpe{p,ip)El{p,ip) 

JG 



2tt 



^1 



dp 



dp^ 



dEn 
dp 



(11) 



where the last integral is taken over the outer surface of 
the cylinder. We choose as unperturbed system for the 
RSE a homogeneous dielectric cylinder in vacuum with 
radius R and refractive index rir (which is neither nor 
1), having 



for 
for 



P^R. 
P> R. 



(12) 



Due to the cylindrical symmetry, the azimuthal index m 
is a good quantum number which takes integer values giv- 
ing the number of field oscillations around the cylinder. 
The unperturbed RS wave functions factorize as 



En{p, = R,n{p, kn)Xyn{'P) , 

where the angular parts are defined by 



(13) 



vr sm{m,tp) if m < , 
Xrniv) = { (27r)-i/2 if m = 0, (14) 

TT^-'^/^ cos(m</j) if m > , 



and are orthonormal according to 



(15) 
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The choice of the wave functions in the form of standing 
waves Eg. (ITU) , instead of the more usual e*""^, is dic- 
tated by the general orthogonality condition defined by 
Eq. without using the complex conjugate. The radial 
components have the form 



Rm{p,k) = A 



J-mirirkp) / Jmin-rkR) 
H„,{kp)/H^{kR) 



for p ^ R , 
for p > R, 



(16) 

in which Jm{z) and Hm{z) = Hm' [z) are, respectively, 
the cylindrical Bessel and Hankel functions of the first 
kind. The wave functions are normalized according to 
Eq. (jlip with the normalization constant 



(17) 



The two boundary conditions at the surface of the cylin- 
der, the continuity of the electric field and its radial 
derivative, produce a secular equation for the RS wave 
number eigenvalues A:„, which has the form 




D^(KR)^Q, 



(18) 



where 



Dm{z) 



■ nrJ'^{nrz)Hm{z) - J^{nrz)H'^{z) , (19) 



and J'm{z) and H'^{z) are the derivatives of Jm{z) and 
Hmiz), respectively. Here z represents a complex argu- 
ment, as opposed to the spatial coordinate used earlier. 

The Hankel function Hm{z) which describes the field 
outside the cylinder and contributes to Eqs. and 
is a multiple-valued function, or in other words is de- 
fined on a Riemann surface having infinite number of 
sheets due to its logarithmic component. However, only 
one of these sheets contains the eigenvalues fc„ , satisfying 
Eq. (jl8l) , which correspond to the outgoing wave bound- 
ary conditions of RSs. This 'physical' sheet of H„i{z) has 
a cut in the complex z-plane along the negative imagi- 
nary half-axis, as shown in Appendix [b1 which in turn 
gives rise to the same cut in the GF. Consequently, the 
Mittag-Leffler theorem used for the spectral representa- 
tion of the GF needs to be modified to include the cut 
contribution, as done in AppendixlBl 

For the TM case treated here, the full GF of the 
homogeneous dielectric cylinder, which is defined via 
Maxwell's equation with a line current source term, is 
given by Gfc = Gk z ® z, in which 

Gfe(p, p') = J2 P'; k)Xrn{^)Xn.{^') , (20) 

m 

and the radial components have the following spectral 
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FIG. 1: (a): Cut poles ka (stars) representing the cut of the 
GF of a homogeneous dielectric cylinder with rir = 2, in the 
complex wave-number plane for m = 0, 11, and 20. Normal 
poles kn (open squares) are also shown, (b): Cut pole density 
o"m(fc) (solid curves) and the cut pole strength </>L'"' (stars), 
for the same values of m. 



representation 



Gm{p,p']k) = ^ 



RmiPi kn)Rm(,P i ^n) 

° Rmjp, k')R^{p',k') 
-ico 2k{k — k') 

RmiPj kn)Rm{p i kn) 

2k(^k k^i) 



<jUk')dk' 
(21) 



derived in Appendix |B] Note that the cut contribution 
to the GF spectrum in the form of the integral in the last 
equation is described in terms of the same functions as 
those used for discrete poles. This implies that the cut of 
the GF can be understood as a continuous distribution 
of additional poles along the negative imaginary half-axis 
with the density 



.ik) 



4(n^ _ l)Jjn{nrkR) 
'K^kD+{kR)D^,{kR) 



(22) 



calculated in AppendixlBl Here D^{z) are the two limit- 
ing values of Dm{z') for z' approaching point z on the cut 
from its different sides Re z' ^ 0. Remarkably, the inte- 
grated density of the cut contribution to the GF is equiv- 
alent to half a normal pole: /^j^Q crm{k)dk — (— 1)™+^/2. 

To numerically treat the cut contribution in the linear 
eigenvalue problem Eq. (|S]), we discretize the integral in 
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Eq. (|2ip into a finite number of cut poles and add cut RSs 
to the basis. These cut poles have non-integer strength 

<j)"a^^ determined by the cut pole density (Jm- The func- 
tion am{k) is purely imaginary and is peaked close to nor- 
mal poles kn as can be seen in Fig.[T]for selected m. In the 
numerical calculations of the present work we have used 
cut pole positions and strengths determined by splitting 
the cut interval [0, — icx)] into ivi™^ regions [9!'"'', ^i+i] 



numbered by a — 1,2. 



(m) 



which are chosen to 



contain an equal weight according to 



\/\a~Jk)\dk = 



N, 



(m) 



VMk)\dk. (23) 



For the numerical results shown later in this section and 
for the chosen values of ivi™\ using the weight v^javij 
in Eq. (|23p was found to give the best accuracy of the 
RSE as compared to other powers of |cr„i|. Each region 
of the cut is represented by a cut pole of the 
GF at fc = fci™'' given by the first moment, 



kam(k)dk 




(24) 



where the cut pole strength 



dm) 



is defined as 



am{k)dk . 



(25) 



An example of cut poles assigned for rn = 0, 11, and 
20 is given in Fig.[IJa). The cut poles contribute to the 
RSE in the same way as the normal poles, and the ma- 
trix elements with the cut RSs are given by the overlap 
integrals Eq. (ID1[) expressed in terms of exactly the same 
functions Eq. (IT6l) as for the normal RSs. In discretiza- 
tion of the linear eigenvalue problem Eq. ^ of the RSE, 
the only modification caused by the cut is that the matrix 
of the perturbation is weighted according to the cut pole 

strengths (j)'^^ , as described by Eq. (|C7p in Appendix [C] 
In the numerical calculation, the total number of poles 
Nt used in Eq. ([S]) determines the computational com- 
plexity of the matrix eigenvalue problem, so that we are 
interested in the number of cut poles in the basis pro- 
ducing the best accuracy for a given Nt- We have in- 
vestigated this numerically for the examples given be- 
low, and found that this is achieved using about 20% cut 
poles in the basis. Only for the homogeneous perturba- 
tion in Sec Hira we used ivi"^ ~ where 7V(") 
is the number of normal poles in the basis for the given 
m, in order to demonstrate the convergence towards the 
exact solution. For all other numerical results we used 



(m) 



0.2iV('^ 



The rest of this section discusses results of the RSE 
for different effective 2D systems. We consider a homo- 
geneous dielectric cylinder of radius R and refractive in- 
dex rir = 2 (e = 4) with different types of perturbations. 



namely a homogeneous perturbation of the whole cylin- 
der in Sec. nil Al a half-cylinder perturbation in Sec. lIIIBl 
a thin-film perturbation in Sec. lIII Gl and a wire perturba- 
tion in Sec. lIII E)] Explicit forms of the matrix elements 
for these perturbations and details of their calculation 
are given in Appendix [Dl 



A. Homogeneous Cylinder Perturbation 

The perturbation we consider in this section is a ho- 
mogeneous change of e over the whole cylinder, given by 



Ae{p,ip) ^Ae0{R-p) 



Ae for p^R, 
for p > i? 



(26) 
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FIG. 2: (a): Perturbed RS wave numbers for the homo- 
geneous perturbation Eq. (|26p calculated via the RSE with 
— 800 (only sine modes are shown). The perturbed poles 
with (-I-) and without (x) the cut contribution are compared 
with the exact solution (open squares). Unperturbed wave 
numbers are also show (open circles with dots). Inset: Di- 
electric constant profile for the unperturbed and perturbed 
systems, (b): Relative error in the calculated perturbed wave 
numbers with (heptagons) and without (triangles) contribu- 
tion of the cut. Relative error for a simulation including the 
cut and improved by extrapolation is also shown (crossed cir- 
cles). 
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with the strength Ae = 4 used in the numerical cal- 
culation. For ((9- independent perturbations, modes with 
different azimuthal number m are decoupled, and so are 
even and odd (cosine and sine) modes given by Eq. (|14p. 
We show only the sine modes here, and use for illustra- 
tion m = 20. The matrix elements of the perturbation 
are calculated analytically and given by Eqs. (jD3l) and 
(ID4p . The homogeneous perturbation does not change 
the symmetry of the system, so that the perturbed modes 
obey the same secular equation Eq. (ITOl) with the refrac- 
tive index rir of the cylinder changed to \/ n1 + Ae, and 
thus the perturbed wave numbers calculated using the 
RSE can be compared with the exact values k\j'^^'^^\ 

We choose the basis of RSs for the RSE in such a way 
that for the given azimuthal number m and the given 
number of normal RSs N we find all normal poles |fc„| < 
knxa^{N) with a suitably chosen maximum wave vector 
kraax{N) and then add the cut poles. We find that as we 
increase N , the relative error \xi, / x^^^^'^^'^ — l| decreases 
as N~^. Following the procedure described in Ref.ls^ 
we can extrapolate the perturbed wave numbers. The 
resulting perturbed wave numbers are shown in Fig. [5] 
The perturbation is strong, creating 3 additional WGMs 
with m = 20 having up to 4 orders of magnitude narrower 
linewidths. For N = 800, the RSE reproduces about 
100 modes to a relative error in the 10~^ range, which 
is decreasing by one or two orders of magnitude after 
extrapolation. The contribution of the cut is significant: 
Ignoring the cut leads to a relative error of the poles in 
the 10~^ range. The fact that the relative error improves 
by 4-5 orders of magnitude after taking into account the 
cut in the form of the cut poles shows the validity of the 
reported analytical treatment of cuts in the RSE, and 
the high accuracy of the discretization method into cut 
poles. 



B. Half-Cylinder Perturbation 

We now consider a bulk perturbation which mixes 
modes with different m. The perturbation is given by 

In our numerical simulation we take Ae = 0.2. The ma- 
trix elements of the perturbation are given by Eqs. (|D5p - 
(ID7p which require a numerical integration. Owing to 
the symmetry of the perturbation, the sine and cosine 
basis modes are still de-coupled, therefore we treat them 
separately, see panels (a) and (b) in Fig. [3] Due to a 
relatively small perturbation (compared to that consid- 
ered in Sec. nil A] ), the mode positions in the spectrum 
do not change much. However, the quality factors Q of 
all WGMs decrease, as the lifetime of the resonances is 
now limited by an additional scattering at the step in the 
dielectric constant of the perturbed cylinder. 



To the best of our knowledge, an analytic solution for 
this perturbation is not available and thus we cannot cal- 
culate the relative error of the RSE result with respect 
to the exact solution. However, we can investigate the 
convergence of the method in order demonstrate how the 
RSE works in this case which is not reducible to an ef- 
fective one-dimensional problem. This is done in Fig. [3] 
showing the perturbed sine and cosine modes for two 
different values of basis size N and in Fig.|4] where abso- 
lute errors M^, are shown for several different values of 
N . Following Ref.[3^ the absolute error is defined here 



as 



maxi=i,2,3 \ ^^^ 



, where are the RS 



wave numbers calculated for basis sizes of Ni « N/2, 
N2 ~ N/V2, N3 « N/^, and N4 = N. The resuhs 
for the cosine and sine modes are quite similar. From 
Fig. m we see that the perturbed resonances are converg- 
ing with increasing basis size. Though, the absolute error 
has some fiuctuations within an order. 

We were able to see the power law of the convergency, 
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FIG. 3: Unperturbed (open circles with dots) and perturbed 
RS wave numbers of (a) cosine and (b) sine modes of a cylin- 
der for a half-cylinder perturbation Eq. (|27p with Ae = 0.2 
and the basis sizes A'^ — 2000 (crosses) and A'^ = 4000 
(hexagons). Only the WGM region is shown. Inset: Diagram 
showing the regions of increased (solid blue) and decreased 
(red and white striped) dielectric constant. 
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the half-cyhnder perturbation Eq. (|27|l as functions of Re>f^, 
calculated via the RSE for different basis sizes A'", for cosine 
(closed shapes) and sine (open shapes) modes. 



in agreement with Ref.[32, and found that the power law 
exponent is approximately —2. We found however that 
owing to the above mentioned fluctuations, the power law 
convergence is not well developed compared to the one- 
dimensional problems considered so far (including the 
example of Sec. lIII A)) . Increasing the basis size N im- 
proves the power-law convergency. This is attributed to 
the larger number of basis states below a given |fcmax-R| in 
effective 2D systems compared to effective ID systems. 
Thus one needs a larger basis in order to approximate 
discrete steps of the basis size by a continuous power 
law. 

To show how a particular perturbed state is created as 
superposition of unperturbed states, we show by a star in 
Fig.[5]one the perturbed WGMs of Fig.[3](a) which, owing 
to the perturbation, increases its linewidth by nearly an 
order of magnitude. The contribution of the basis states 
is visualized by circles of a radius proportional to |c„y |^/'^, 
which are centered at the positions of the wave vectors 
kn in the complex /c-plane. The expansion coefficients 
Cnu decrease quickly with the distance to the spectral 
position of the perturbed mode >c^, with the dominant 
contribution coming from the nearest unperturbed RS, a 
typical feature of perturbation theory in closed systems. 
Importantly, this demonstrates that if we are interested 
in the modes within a small spectral region, we can limit 
the basis in the RSE to states close to that region. This 
result is crucial for the application of the RSE to effec- 
tive 3D systems which have even larger numbers of basis 
states below a given l/cmax^l, as one can significantly re- 
duce the number of basis states needed to calculate the 
perturbation of a mode of interest to a given accuracy. 
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FIG. 5: Contributions of the basis RSs (black circles) to 
a given perturbed RS (blue star), calculated for the co- 
sine modes and the half-cylinder perturbation Eq. lf27)l , using 
N = 4000. All circles and the star are centered at the posi- 
tions of the corresponding RS wave numbers in the complex 
fc-plane. The radius of the black circles is proportional to 
|c„^|^'''^. A key showing the relationship between circle radius 
and |c„^p is given by the red circles. 



C. Thin-Film Perturbation 

We now move from the bulk perturbations towards the 
case of a thin film embedded in the cylinder, corresponds 
to a line perturbation in the effective 2D system. The 
perturbation we consider in this section is given by 



Ae(p, V?) = hAs 



0{R-P) 



P 



(28) 



see the inset in Fig. [51 In our numerical simulation we 
take the strength of the perturbation hAe = —O.IR. 
Physically this perturbation corresponds to a thin metal 
film of uniform negative dielectric constant + Ae and 
width h much narrower than the shortest wavelength in 
the basis. The perturbation leaves the sine modes of the 
unperturbed cylinder unchanged. Hence we only include 
cosine modes into the basis. The perturbation matrix 
elements are given by Eq. (jD8|) . 

To our knowledge an analytic solution for this per- 
turbation is not known and we therefore calculate the 
absolute error as in Sec. lIII Bl Fig. [5] shows the resulting 
RS wave numbers and absolute errors for this thin-film 
perturbation. We see in Fig.|6ljb) that the convergence 
of the RSE is slower than in the case of the half-moon 
perturbation. This is expected as the thin film has no 
geometrical effect on the wave- vector ky, giving higher 
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FIG. 6: (a): Unperturbed and perturbed RS wave numbers 
of cosine modes for a thin-film perturbation given by Eq. (|28p 
with hAe — —O.IR, calculated via the RSE with the basis 
sizes N = 2000 (crosses) and = 4000 (hexagons). The un- 
perturbed RSs are shown as open circles with dots, (b): Ab- 
solute errors as functions of Re x-^ calculated for different 
basis sizes A'^ as labeled. Inset: Sketch showing the location 
of the thin metal film perturbation as a red line inside the 
unperturbed cylinder. 

contributions of basis states with large ky , similar to the 
results in ID with a delta scatterer perturbation reported 
earheii^ We have found that the power law exponent in 
this case is approximately —1. 



D. Thin- Wire Perturbation 

As last example we consider a dielectric cylinder per- 
turbed by a thin-wire perturbation which is represented 
by small disk of radius b centered at the point d on the 
X-axis {if = 0). We do not use here a delta perturbation, 
in order to compare it with an analytic solution available 
in the literaturei^ The perturbation is defined as 

As{p) = Ae9{b-\p~d\), (29) 
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FIG. 7: (a): Unperturbed and perturbed RS wave numbers 
of cosine modes for a thin-wire perturbation given by Eq. (|29p 
with d = 0.8R, b = O.OOli?, and Ae = 100, calculated via the 
RSE with the basis size N — 2000 (crosses) and compared 
with the analytic approximation of Ref.lsJ (empty squares). 
The unperturbed basis states are shown as open circles with 
dots, (b): Absolute errors in RSE as functions of Kexi, 
calculated for different basis size A'^ and the absolute difference 
between the RSE and the analytic approximation (crosses x). 
Inset: Sketch showing the location of the wire as a red dot 
inside the unperturbed cylinder. 



and we choose d = |d| = 0.8R, b = O.OOli?, and Ae = 100 
(the unperturbed system is the same as before having 
Ur — 2.0). This perturbation leaves the sine and cosine 
modes decoupled, and the sine modes approximately un- 
changed (strictly for d — 0). Therefore we show here 
the perturbation of the cosine modes. The RSE pertur- 
bation matrix elements are given by Eq. (IDIOI) . The re- 
sulting RS wave numbers are shown in Fig.[7](a) together 
with the analytic approximation, demonstrating a good 
agreement. 

The absolute errors are shown in Fig.[7](b). We 
see that the convergence in the case of a thin wire is 
even slower than for the thin-film perturbation shown 
in Sec. nil Cl This is expected as the thin wire has no 
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geometrical effect on both and ky, giving higher con- 
tributions of the basis states with large |fc|. We found 
that within the basis sizes investigated, the power law 
is not well developed, but for weaker (smaller |Ae|) or 
more spatially extended (larger b) perturbations a better 
convergence is observed, as expected. 

The analytic solution of Ref.ls^ for a point-like scat- 
terer in a 2D disk is not strict in any physical system. 
In the case of a delta scatterer, the secular equation 
is logarithmically divergent and thus cannot be used, 
while the accuracy of the model for a finite size scat- 
terer relies on a number of approximation a'^^'^'^ which 
require \nrXuh\ <C 1, + Ae| ^ 1, and also 

|Re>tfi,| ^ |Im>fi,|, i.e. having a large Q. In addition to 
this, the point-like perturbation should not be too close 
to the edge of the disk, i.e. \nrX,y{R — d)\ <^ 1. While 
we have chosen our parameters to be suitable for these 
approximations, we do not have a quantitative estimate 
of the error. Nevertheless, the comparison in Fig.[7](a) 
of the RSE calculation with the analytic solution demon- 
strates a good agreement which is improving as we move 
closer to the origin in the complex /c-plane, as detailed in 
Fig.[7](b) where the absolute difference between the two 
calculation is shown. 



IV. SUMMARY 

We have applied the resonant state expansion (RSE) 
to effective two-dimensional (2D) open optical systems, 
such as dielectric micro-cylinders and micro-disks with 
perturbations. We have found and treated a cut of the 
Green's functions (GFs) of effective 2D systems - a fea- 
ture which to our knowledge has not been mentioned in 
the literature but turned out to be crucial for the RSE as 
the states on the cut contribute to the completeness of 
the basis of RSs needed for the accuracy of the method. 
We have detailed the formulation of the RSE for a gen- 
eral 3D case taking into account the vectorial nature of 
the electro-magnetic field and tensor form of the GF and 
shown in detail how the theory is applied to effective 2D 
systems for which states on the cut are introduced and 
discretized for the numerics. 

Using the analytically known basis of resonant states 
(RSs) of an ideal homogeneous dielectric cylinder - a 
complete set of eigenmodes satisfying outgoing wave 
boundary conditions - we have treated different types of 
perturbations, such as half- cylinder, thin-film and thin- 
wire perturbations. For all of these perturbations, the 
perturbed systems are not reducible to effective ID ones, 
so that the present work demonstrates the applicability 
of the RSE to general effective 2D perturbations which 
mix all basis modes. We investigated the convergency 
for these perturbations and compared the RSE results, 
where it was possible, with available analytic solutions. 
In particular, we have made such a comparison for a ho- 
mogeneous perturbation of a cylinder, which is reducible 
to an effective ID system, and for point-like perturba- 



tion of a disk which presents an essentially 2D system 
with mixing of all kind of modes in the given polariza- 
tion of light. In both situation we have found agreement 
between the RSE and the known analytic solutions. 
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Appendix A: Spectral representation of the Green's 
function of an open system 

The Green's function (GF) of an open electromagnetic 
system is a tensor Gfc which satisfies the outgoing wave 
boundary conditions and the Maxwell wave equation ([T]) 
with a delta function source term, 

-VxVxGfc(r,r') + /c^£(r)Gfc(r,r') = i(5(r-r'), (Al) 

where 1 is the unit tensor and k = cj/c is the wave vector 
of the electro-magnetic field in vacuum determined by the 
frequency w, which can be real or complex. Physically, 
the GF describes the response of the system to a point 
current with frequency w, i.e. an oscillating dipole. Using 
the reciprocity theorem?^ the relation 

diGfe(ri,r2)d2 = d2Gfc(r2, ri)di (A2) 

holds for any two dipoles di_2 at points ri^2 oscillating 
with the same frequency. Therefore Gfe(r,r') is a sym- 
metric tensor. Assuming a simple-pole structure of the 
GF with poles at k = qn and taking into account its large- 
k vanishing asymptotics, the Mittag-Leffier theorem^i^ 
allows us to express the GF as 

G..(r,r') = ^%(^. (A3) 

Substituting this expression into Eq. (|A1[) and convolut- 
ing with an arbitrary finite field D(r) over a finite volume 
V we obtain 

^ -Vx VxF„(r) + p£(r)F„(r) _ 
^ k-qn 

n 

where F„(r) = Q„(r, r')D(r')dr'. Taking the limit 
k ^ qn yields 

-V X V X F„(r) + g2e(r)F„(r) = . 

Due to the convolution with the GF, F„(r) satisfies the 
same outgoing wave boundary conditions Eq. Then, 
according to Eq. ([!]), F„(r) oc E„(r) and gn = fc„. Note 
that the convolution of the kernel Qn(r, r') with different 
functions D(r) can be proportional to one and the same 
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function E„(r) only if the kernel has the form of a direct 
product: 



Q„(r, r') = E„(r) ® E„(r')/u;„ , 



(A4) 



in which the vectorial nature of the field and the symme- 
try of the kernel Eq. (IA2[) were taken into account and 
a normalization constant w„ was introduced. The direct 
vector product (g) is defined as c(a(g) b)d = (c • a)(b • d), 
for any vectors a, b, c, and d. 

The asymptotics of the GF for /c — > oo following from 
Eq. (jAll) requires that Gk oc k~^, which for the spectral 
representation Eq. ()A3|) provides the sum rule)^ 



(A5) 



5^Q4r,r')=0. 



Using Eq. and Eq. (|X5|) . together with Eq. (jM 

yields 



Gfc(r,r') 



E- 

n 



,(r)(8)E„(rO kn 



(A6) 



The normalization constants Wn can be determined from 
the normalization condition Eq. (|4]) . We have shown for 
specific systems with analytic solutions^ that 



'Zkji 



(A7) 



These systems were ID or 3D with spherical symmetry, 
where the normalization constants Wn have the meaning 
of the Wronskian derivatives at k = kn- The deriva- 
tion for 2D systems with cylindrical symmetry is similar. 
However a specific property of 2D systems is the presence 
of a continuum in the RS spectrum, which is a cut of the 
GF in the complex A:-plane. The spectral representation 
of the GF includes this cut contribution as integral: 



,^ _ W; E„(r)0E„(rO 
2fc(fc-fc„) 



(A8) 



The cut is discussed in more detail in Sec. lIIII and in Ap- 
pendix |b1 where we also show the validity of Eq. (|A7|) for 
a homogeneous cylinder. 

Substituting the GF of Eq. ([Ml) into Eq. and us- 
ing the sum rule Eq. (|A5[) leads to the closure relation 

^^E„(r)®E„(r') = i5(r-r') (A9) 



which expresses the completeness of the RSs, so that any 
function can be written as a superposition of RSs includ- 
ing the contribution of the cut. Due to the surface term 
in the orthogonality condition Eq. ([3]) , RSs form an over- 
complete basis,— i.e. any RS wave function E„(r) can 
be written as a superposition of the wave functions of 
other RSs of the basis. 



Appendix B: Green's function of a homogeneous 
cylinder 

The TM component of the GF of a homogeneous cylin- 
der in vacuum satisfies the following equation 



Vl + e{p)k^ Gk{p,p') = S{p-p') (Bl) 



with 



Using the angular basis Eq. the GF can be written 
as 

Gkip, p') = -i= E ^™('°' P'-' k)Xm{v)Xm{v') , (B3) 



similar to Eq. ([20l) . Note that we redefined here the radial 
part as Gm{p, p'jk) — ^/ pp'Gm{p, p'; k) which satisfies 



dp^ 



- 1/4 



+ Pe{p) 



GUp,p'-k) = 5{p-p'). 



(B4) 

Using two linearly independent solutions fm{p) a-nd 
9m{p) of the corresponding homogeneous equation which 
satisfy the asymptotic boundary conditions 

/„(p)(xp"+i/2 for p^O, 
gm{p) oc e'^'^P for p — > oo . 



the GF can be expressed as 
Gm{p,p'; k) 



Ini{p<)gni{p>) 



(B5) 



in which p< = min{/3, p'}, p> ~ max{/3, p'}, and the 
Wronskian W{f, g) — fg' — f'g . For TM polarization, a 
suitable pair of solutions is given by 



where 



= V~p- 



Jmirirp) , 



p^R, 



am-Jmip) +bmHmip) , P>R, 
CmJminrP) + a„iH„i{nrp) , p ^ i? , 



Hmip) , 



P> R, 



am{k) = [nrJ^{nrx)Hmix) - Jm{nrx)H!^{x)]Trix/2, 

bm{k) = [j!^^{x)Jrn{nrX) - nrJm{x)J^{nrX)]Trix/2 , 

Cm{k) = [H',^{x)Hjn{nrx) ~ nrHrn{x)H'^{nrx)]T:ix/2 
with X ~ kR. The Wronskian is calculated to be 

W{frn,gm) ^ 2ia„i{k) / TT = -xD,n{x) , 

with Dm{x) defined in Eq. (fT9| . Inside the cylinder, the 
GF then takes the form 



G„i{p,p';k) = Tp\/pp' 

Cm{k) 



am{k) 



J„i{nr kp^ ) Hrn {rirkpy ) 



Jm{nr kp^ )Jm{nrkpy) 



(B6) 



The GF has simple poles fc„ in the complex fc-plane whic 
are the wave vectors of RSs, given by am(fc„) = 0, a 
equation equivalent to Eq. ([TSl) . The residues Res„ ( 
the GF at these poles are calculated using 



jjja„i(fc) 



2z /t j2 



7r(n2 - l)[knRJrn{nrkr 



(B^ 

In addition to the poles, the GF has a cut in the con 
plex fc-plane along the negative imaginary half-axis. Tf 
cut is due to the Hankel function Hm{z) which describe 
the field outside the cylinder and contributes to Eqs. ([T6l 
([TSl) and Eq. (jB6p . and is not uniquely defined. Indee( 
it can be expressed as^ 



Hm{z) = Jm{z) + iN„i{z) , 

using a multiple- valued Neumann function 



(B^ 



N^{z) = N^{z) + -J^{z)\n^, (EC 

where iVm(z) = z™Fm{z'^) is a single-value 
polynomial^ while In z is a multiple- valued funi 
tion defined on an infinite number of Riemann sheet 
We have verified that only one such sheet provides tf 
asymptotics Hm{z) oc exp(i2)/Y^ for z — > oo, which 
required for the RS wave functions outside the cylindc 
to satisfy the outgoing wave boundary conditions Eq. ([2] 
This 'physical' sheet has a cut going from the branc 
point at z = to infinity, and the position of th 
cut is not arbitrary. To find the cut position let us„ 
the symmetry of the RS wave numbers, k-n = —k^, 
discussed in SeclHl Let us also, using properties of 
cylindrical functions^i^ bring the secular equation (fTSl) 
to the form 

J„i+linrZ)H„i-liz) = J,n-i{nrz)H,n+i{z) , (BIO) 

in which z — knR. We note that if z = /c„i? is a complex 
solution of Eq. (jBlOp . then —z* is also a solution of the 
same equation. We take two equations, one is the con- 
jugate of Eq. (jBlQ]) and the other is Eq. (IbTO]) itself but 
taken with the argument — z*, and add them up. Sub- 
stituting there Eqs. (|B8|) and (jB9p and using the facts 
thati^ 



[Jm{z)\* = Jm{z*) = 
[7V„(z)]* = Nm{z*) 

we arrive at the condition 



(-l)"J„,(-z*), 
= (-l)'"iV„(-z*), 



ln(- 



(Inz) 



(Bll) 



which is fulfilled, for any z, only if In z [and consequently 
Hm{z)\ has a cut along the negative imaginary half-axis. 

Owing to the cut of the Hankel function Hm{z) the 
GF also has a cut along the negative imaginary half- 
axis in the complex fc-plane, so that on both sides of the 




>Re(k') 



FIG. 8: Sketch showing the contour of integration in Eq. (|B14|) 
as well as poles (black dots) and the cut (blue dashed line) of 
the GF in the complex fc'-plane. An extra pole at fc' = fc is 
shown by a red dot. 



cut Gm takes different values: one the right-hand 
side and G" one the left-hand side of the cut. The step 
AGm = G^ — G~ over the cut can be calculated using 
the corresponding difference in the Hankel function: 

AH„,{z) = H+{z) - H^{z) = 4J^(z) . 

The result is 

~ TT / 

AG„i(p, p';fc) = —y^ pp' Jrn{nrkp^)Jrn{nrkpy)AQjrt{k) 



where 



AQ™(fc) 



(B12) 



nkRj D+{k)D:;^,{k) 
(B13) 

with Dmik) given by Eq. ([T9| . 

Let us now use the residue theorem for the function 
G,n{p,p'',k')/{k — k') integrating it in the complex fc'- 
plane along a closed contour consisting of three parts, 
see Fig.lHl A large counter-clockwise circumference with 
a radius tending to infinity, two straight lines circum- 
venting the cut and approaching it from both sides, and 
a small clockwise circumference around the origin with 
a radius tending to zero. Since the GF behaves as fc~^ 
at large values of fc and takes finite values or logarith- 
mically diverges at fc = 0, both large- and small-circle 
integrals vanish, so that the only remaining integrals are 
those which are taken along the cut: 



GUp,p';k') 



2TTi 



k-k' 
■sr-^ Res, 



dk' ^ 



G+dk' 



Q k ' k^ 
2TTiG„i{p,p';k) . 



k-k' 
(B14) 
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Note that in second part of the above equation we have 
made use of the residue theorem, expressing the closed- 
loop integral in the left-hand side in terms of a sum over 
residues at all poles inside the contour. Using Eq. (|B14[) 
the GF can be expressed as 



where 



Gm{p,p'; k) 



Res„ 



^ k 



1 

"27ri 



l^Gm{p.p'\k')dk' 
k-k' 



(B15) 

which is a generalization of the Mittag-Lefher theorem. 
This is used in Sec. lIIII when applying the RSE to 2D 
systems with a cut of their GF. The residues Res„ of the 
GF contributing to Eq. (jB15[) are calculated as 



Res„ 



2i 



^ Jni{nrkp<)Jm{nrkpy)rra[kn) , (B16) 



with rm{kn) found in Eg. (jB7p . Given that the spa- 
tial dependence of the GF, as described by Eqs. (jBlsp . 
()B16p . and (jB12p . is represented by products of the RS 
wave functions Rm{p,kn) and their analytic continua- 
tions Rmip, k) with fc- values taken on the cut, we arrive 
at the GF in the form of Eqs. .(^1]) and ([^ which are 
then used in the RSE. 



Appendix C: RSE with a cut 

This section provides some details on how the RSE can 
be used in practice when the GF has a cut, and in par- 
ticular, how the cut discretization, producing cut poles, 
modifies the linear matrix eigenvalue problem Eq. ([8|) , the 
central equation of the RSE. 

To make our consideration as general as possible we 
use the 3D version of the Dyson equation ([5]) as a start- 
ing point and substitute into it spectral representations 
Eq. (jA6l) of the unperturbed and perturbed GFs. Equat- 
ing the residues at the perturbed poles (i.e. integrating 
along the infinitesimal circumference around each pole) 
we obtain^ 



E„(r)/E„(rO-£.(rOA£(rOdr' 

2{kn/>c^ ~ 1) 



(CI) 



where both pole and cut contributions are included. The 
cut has a continuous contribution. Therefore, in numeri- 
cal calculation, we discretize the cut representing it by a 
finite number of cut poles chosen in an optimum way as 
described in Sec. lIIII For an arbitrary function F{k), 

^i^„ = J2 ^(^") + r F{k)a{k)dk « P^-^'t'^ 

(C2) 

where the combined index n is used to denote both real 
poles kn and cut poles ka simultaneously. The weighting 
factors (j)n come from the weight function a{k) introduced 
in Eq. ([22l) and are defined as follows 



1 for real poles, 
for cut poles, 



(C3) 



9c + l 



a(k)dk 



(C4) 



is an integral of the weight function over an interval 
[<ion<la+i\ within which a cut pole ka is chosen. The 
method of choosing this interval and the positions of the 
cut poles specific to this work are described in Sec. lIIII 
by Eqs. (|23|) and (|24|) . Note that we have dropped here 
the azimuthal index m for simplicity, however the dis- 
cretization of the cut is generally different for different 
m leading to m-dependent cut poles and their weighting 
factors (pn- 

Now, discretizing Eq. (jCll) in accordance with 
Eq. (jC2l) , substituting the expansion 



(C5) 



= ^6„^E„(r) w ^0fi6^^Efi(r) 



of the perturbed RS wave function into it, and equating 
coefficients at the same basis functions, we obtain 



1 



2{kn/>t, - 1) ^ 

n' 

Introducing new expansion coefficients 



(C6) 



— h I ^" 

the matrix eigenvalue problem takes the form 

Vnn' I (f>n4>n' \ _ _ 1 




(C7) 



which is a discretized version of Eq. . The matrix el- 
ements Vnn' are defined by Eq. ([9]) with n and n' num- 
bering both normal and cut poles. They are calculated 
below in Appcndix[D] for various types of perturbation. 



Appendix D: Matrix elements for various 
perturbations in 2D 

In this section we give explicit expressions for the ma- 
trix elements Vnn' of the specific perturbations consid- 
ered in this paper. As a starting point we use the follow- 
ing general formula for the matrix elements of an arbi- 
trary perturbation Ae(p, ip) inside the cylinder of radius 
R: 



Vnn' 



2ir 



XmW)X^n'W)dip 



X / Ae{p,(p)Rm{p,kn)Rjn'{p,kn,)pdp , (Dl) 

^0 

in which Rm and Xm are the eigenfunctions of the homo- 
geneous cylinder given by Eqs. (ini) - pT|) . 
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1. Homogeneous cylinder perturbation 



The homogeneous perturbation Eq. ((26)) does not mix 
different m- values. The matrix elements between RS 
with the same azimuthal number m are given by the ra- 
dial overlap integrals 



(D2) 



Vfin' = Ae / Rm{p,kn)Rm{p,kn')pdp 

Jo 

yielding for identical basis states {n — n') 
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and for different basis states (n ^ n') 
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2. Half-cylinder perturbation 

The most efficient way of calculating the matrix el- 
ements of the perturbation Eq. (|27|) is to calculate the 
angular parts of the integrals analytically and the radial 
parts numerically. The matrix elements have the form 



(D5) 



in which the angular overlap integrals Pmm' are vanish- 
ing when taken between modes of different parity, i.e. be- 
tween sine and cosine modes, see Eq. ([T4l) . and between 
same parity modes corresponding to azimuthal numbers 
m and m' of different parity. The non- vanishing integrals 
are given by 

/7r/2 r.Z-K/2 
X7n{v)X7n'{v)d^ - / X7n{'f)X77i'{v)dV 
-Txjl J-k/2 

= S,„S„i'(-0m-m' ± V'm+m') (D6) 

with -|- (— ) corresponding to cosine (sine) modes and Sm 
and iprn defined as 

TT^^/-^ for m 7^ , 
(27r)-i/2 for m = 0, 

^™ = [i-(-ir]^^^i^^. 



The radial part of the matrix elements of the perturba- 
tion is given by the integrals 
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which we calculate numerically. 



3. Thin-film perturbation 

The matrix elements of the perturbation Eq. (|28p are 
jiven by the integrals 



Vnn'=hAexlM R7n{pM)R7n'{pM')dp, (D8) 
Jo 



similar to Eq. (jD7p , which are calculated numerically. 



4. Thin-wire perturbation 



The RSE perturbation matrix elements for this system 
are calculated by summing / same-strength delta scatter- 
ers on a square grid covering a circle. The perturbation 
Eq. ([29]) is thus modeled by 



Trfe^ 1 

Ae w Ae— V -5{p ~ p7)5{if - ip^) . 
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The matrix elements then have the form 



Vnw = A£^^Ef,ip^,ip^)Ef,'{p7,ip7) (DIO) 



with En{p, ip) given by Eq. P^ . 
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